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We study the behavior of hybrid stars using an extended hadronic and quark SU(3) non-hnear 
sigma model. The degrees of freedom change naturaUy, in this model, from hadrons to quarks as the 
density /temperature increases. At zero temperature, we reproduce massive neutron stars containing 
a core of hybrid matter of 2 km for the non-rotating case and 1.18 km and 0.87 km, in the equatorial 
and polar directions respectively, for stars rotating at the Kepler frequency (physical cases lie in 
between). The cooling of such stars is also analyzed. 
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I. INTRODUCTION 

As the densest bodies in the universe, neutron stars 
are one of the best suited candidates to look for quark 
matter. Studying when and how the deconfinenient from 
hadronic to quark matter occurs is crucial for the under- 
standing and identification of such phenomena. Usual 
approaches for hybrid neutron stars consist of two differ- 
ent models with separate equations of state for hadronic 
and quark phases (see e.g. |1|), connected at the point for 
which the pressure of the quark phase exceeds that of the 
hadronic one. Within our approach we employ a single 
model for the hadronic and quark phases, avoiding the 
need for two separate equations of state. This approach, 
that presents a more natural transition from hadronic to 
quark matter, allows us to follow the spin evolution of a 
hybrid star with a single equation of state. In this pa- 
per we will investigate the structural changes that follows 
the stellar spin-down and the subsequent effects on the 
cooling of the object. 

The SU(3) non-linear sigma model introduces baryons 
and quarks as flavor-SU(3) multiplets. Baryons and 
quarks obtain their masses through their coupling to the 
scalar fields of the theory (with an additional coupling 
to the Polyakov loop as discussed below) via sponta- 
neous symmetry breaking. We include the quark de- 
grees of freedom in the hadronic model in analogy to 
the PNJL model [| by using an effective field that 
can be related to the QCD Polyakov loop, defined via 
<i> = iTr[exp [i J dTA4)], where A4 = iAg is the tempo- 
ral component of the SU(3) gauge field. The effect of the 
field is to suppress quarks in the low density/temperature 
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regime and baryons at high densities and temperatures, 
respectively. 

This paper is divided as follows: in section II we review 
the properties of the model used for the composition and 
equation of state; in section III we present our results 
which encompass the structure of rotating and spheri- 
cally symmetric compact stars, and cooling effects; and 
in section IV our conclusions are presented. 



II. THE MODEL 

The Lagrangian density for the sigma-type model, in 
mean field approximation, is given by 



L — LKin + Lint + Lself + LsB — U, 



(1) 



where besides the kinetic energy term for hadrons, 
quarks, and leptons, the terms 



Lint = - Li '<Pi[lo{9iici^ + gi4><l> + 9ipT?.p) + M*]iIj„ 

Lself = -i(m2 _|_ ^2^2 ^2^2) 

+ ko{a^ + e + S^) + hia^ + e + S^f 
+ + f + 3(72,52 + C*) + ksia^ - 6^)C 



(2) 



-|- In 



LsB = mlf-^a + [V2mlfk - -^mlf^) (, 



(3) 



(4) 



represent the baryon (and quark) - meson interactions, 
meson self-interactions, and an explicit chiral symme- 
try breaking term that is responsible for producing the 



2 



'1 


1 


— p,n 
.... ^ 

- u,d 
■ ■ s 





























,1 


1 , 1 



1000 1200 1400 

|Xg (MeV) 



1600 



2.5 




1.5 



0.5 



J0751+1807, M = 1.26 + 0.14 



2 ^ 

hadronic star ' ^ ^ 

hybrid .star (local c. n.) i gS 



hybrid star (global c. n.) 



12.8 13 13.2 13.4 



10 



11 12 
R(km) 



13 14 



FIG. 1. (Color online) Effective normalized mass of different 
species as a function of baryonic chemical potential for star 
matter at zero temperature Q. 

masses of the pseudo-scalar mesons. We will discuss the 
potential U further along this paper. The model has 
a SU(3) flavor symmetry, and the index i denotes the 
baryon octet and the three light quarks. In our calcu- 
lations we take into account the following mesons: the 
vector-isoscalars uj and 0, the vector-isovector p, the 
scalar-isoscalars a and C (non-strange and strange quark- 
antiquark states, respectively), and the scalar-isovector 
5. The coupling constants of the model can be found 
in reference They were fitted to reproduce the vac- 
uum masses of the baryons and mesons, nuclear satu- 
ration properties (density po = 0.15 fm~^, binding en- 
ergy per nucleon B/A = —16.00 MeV, nucleon effective 
mass M'^ = 0.67 Mn, compressibility K = 297.32 MeV), 
asymmetry energy {Esym — 32.50 MeV), and reasonable 
values for the hyperon potentials {U\ = —28.00 MeV, 
Us = 5.35 MeV, Us = -18.36 MeV). The vacuum expec- 
tation values of the scalar mesons are constrained by re- 
producing the pion and kaon decay constants. A detailed 
discussion of the purely hadronic part of the Lagrangian 
can be found in [l-Q. 

The effective masses of the baryons and quarks are 
given by 

^Ib = 9bc,<7 + gssTaS + gscC + Mq^ + 35$$^, (5) 
M* = gqao- + gqSTsS + gqcC + + .gg$(l - $), 

(6) 

where AIq is equal to 150 MeV for nucleons, 354 MeV for 
hyperons, 5 MeV for up and down quarks and 150 MeV 
for strange quarks. 

Equations ([5]) and ([6]) show that as the field $ increases 
(with the increase in density/temperature) baryons are 
supressed, giving way to the quark phase, effectively 
modeling the QCD deconfinement phase transition. The 



FIG. 2. (Color online) Mass-radius diagram for the model 
investigated in this paper. The horizontal lines represent ob- 
served pulsar masses ([^-[lEl and references therein) . The 
band delimited by the two thick horizontal lines (orange lines 
online) represent the range of pulsars masses observed, ac- 
counting for the error of the highest and lowest observed mass. 



opposite is true for low values of <i> (at low den- 
sity / temperature) . 

The effective normalized masses of baryons and quarks 
are shown in Fig. [1] Since the coupling constants in the 
<i> term of the effective mass formulas are high but still 
finite, the effective masses of the species not present in 
each phases are large but also finite. 

In analogy to the PNJL model, we define the potential 
U for <i> as 



+ a3ro4log(l-6$2 + 8$3_3$4)_ 



In our case, U{^) is a simplified version of the potential 
used in 0, § and adapted to include terms that depend 
on the chemical potential. These two extra terms are not 
unique, but the most simple natural choice of extending 
the potential. The corresponding parameters are chosen 
to reproduce the main features of the phase diagram at 
finite densities. The coupling constants for the quarks 
can be found in reference [3,], and are chosen to reproduce 
lattice data as well as known information about the phase 
diagram. The lattice data includes a first order phase 
transition at T = 270 MeV, and a pressure function P(T) 
similar to refs 0, [1| at p = for pure gauge (for the 
quenched case without hadrons and quarks). 
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FIG. 3. (Color online) Population (baryonic density for dif- 
ferent species as a function of baryonic cliemical potential) for 
star matter at zero temperature using global charge neutrality 
i- 



FIG. 4. (Color online) Equation of state (pressure as a func- 
tion of energy density) for star matter at zero temperature 
using local and global charge neutrality ^J. 



III. RESULTS 
A. Deconfinement to Quark Matter 

In our model the quarks are suppressed in the hadronic 
phase and the hadrons are suppressed in the quark phase, 
up to = 1700 MeV for T = 0. This behavior is due to 
the fact that the coupling constants in the $ term of the 
effective mass formulas are high but still finite, so in prin- 
ciple at very high chemical potential the threshold can be 
reached a second time for hadrons. This threshold, which 
is higher than the density in the center of neutron stars, 
establishes a limit for the applicability of the model. The 
hyperons, despite being included in the calculation, are 
suppressed by the appearance of the quark phase. Only 
a very small amount of Lambdas appears immediately 
before the phase transition. The strange quarks appear 
after the other quarks with relatively low abundance. 
The density of electrons and muons is significant in the 
hadronic phase but not in the quark phase, since the 
down and strange quarks are negatively charged, reduc- 
ing the need for leptons to maintain charge neutrality. 

Focusing on the high-density/low-temperature region 
of the phase diagram we are able to calculate neutron 
star masses and radii solving the Tolmann-Oppenheimer- 
Volkof equations [l^, [l3|. The solutions for hadronic 
(same model but without quarks) and hybrid stars are 
shown in Fig. [51 where besides our equation of state for 
the core, a separate equation of state was used for the 
crust [18|. The horizontal lines in Fig. [5] represent ob- 
served masses of some prominent millisecond pulsars (0- 
[Tsl and references therein) . To avoid cluttering the graph 
we only indicated the name of the pulsars that establishes 
the lowest and highest observed masses: J0751+1807 



with mass M/Afsun = 1-26 ± 0.14 [Tl| and J0437-4715 
with M/Msun = 1.76±0.2 [i,[i2| respectively. Any plau- 
sible equation of state must be able to produce neutron 
stars within the range delimited by these two objects, as 
is the case for the model investigated in this paper. Our 
model predicts a maximum neutron star mass of 2.1 Mq 
considering local charge neutrality, which as shown in 
Fig. [21 is in agreement with the observational constraints 
shown. Comparing to the hadronic star sequence, the hy- 
brid sequence has a much sharper peak. This peak (which 
denotes the most massive stable star in the sequence) sig- 
nals the phase transition into quark matter at the core of 
the star. Because the equation of state for quark matter 
is much softer than the one for hadronic matter, the star 
becomes unstable at the point where the central density 
is higher than the phase transition threshold. 

There is still another possible option for the configu- 
ration of the particles in the neutron star ^Ift] . If instead 
of local we consider global charge neutrality, we find a 
mixture of phases even for zero temperature as discussed 
in [20]. This possibility, which is a more realistic ap- 
proach, changes the particle densities in the coexistence 
region making them appear and vanish in a smoother 
way (Fig. [3]). Therefore, the maximum mass allowed for 
the star is slightly lower in this case than in the previ- 
ous one, as can be seen from the dotted line in Fig. [51 
however, this possibility allows stable hybrid stars with 
a small amount of quarks. The mixed phase constitutes 
the inner core of the star up to a radius of approximately 
2 km. The equation of state for both cases is shown 
in Fig. [4l The large jump in the pressure for the local 
charge neutrality case explains why the neutron stars be- 
come immediately unstable after the phase transition in 
this configuration. 



4 




200 300 400 500 600 700 

(MeV/fm^) 



FIG. 5. (Color online) Gravitational mass as a function of cen- 
tral density. The lowermost sequences represent static stars 
(n = Hz) , whereas the higher most represents stars rotating 
at their Kepler frequency. The curves connecting the Kepler 
frequency to the static sequence are for stars with constant 
baryonic mass. These sequences represent the evolution of 
Kepler frequency stars to non-rotating objects. 



B. Rotational Effects 

The results shown in Fig. [2] are for static neutron stars 
without including rotational effects. The rotational na- 
ture of pulsars (which in some cases may be rotating with 
frequencies as high as 700 Hz) warrants the investigation 
of the structure of rapidly rotating neutron stars, which 
is considerably more cornplicated than that of static ob- 
jects. In References d, [2l| the rotational effects of the 
hadronic part of our model were investigated by means 
of the improved Hartle-Thorne pertubative method [25|. 
Here we will extend this research by performing exact 
rotational calculations and by including the quark-phase 
described in previous sections. The numerical method 
used for the solution of Einstein's field equations, and for 
the stellar structure of rapidly rotating neutron stars is 
based on the KEH method [2j|, which basically consists 
of expanding the metric functions in terms of Green's 
functions, which can be iteratively integrated, allowing 
us to calculate the structure of the star. This method 
has been expanded by several authors, and details can 
be found in References (23 . [25| . 

In Fig. [S] we show the mass as a function of central 
density. The lowermost sequence represents static stars 
(17 = Hz), whereas the highest one represents stars ro- 
tating at their Kepler frequency. The Kepler frequency, 
or mass shedding frequency, is the maximum frequency at 
which a compact star may rotate. When rotating above 
this frequency an object would shed mass at its equa- 
tor. This quantity therefore sets an absolute limit for 
the frequency of compact stars. The curves connecting 
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FIG. 6. (Color online) Gravitational mass as a function of 
frequency for the stars of Fig. [5] The stars from the static 
sequence are on the y-axis of the graph. 

the Kepler frequency to the static sequence indicate stars 
with constant baryonic mass. These sequences represent 
the evolution of Kepler frequency stars to non-rotating 
objects. This is better illustrated in Fig. [BJ which shows 
the gravitational mass as a function of frequency for the 
constant baryon mass stars. Fig. [B] also shows that the 
gravitational mass may be as high as 3% higher for a 
star rotating at the Kepler frequency (compared to a 
non- rotating star with same baryonic mass). The cen- 
tral density, on the other hand, may be decreased by 
as much as 38%, as shown in FigjSl It is important to 
stress that the constant baryon mass sequences shown in 
Figs. [S] and [B] represents the spin-down evolution of an 
isolated compact star (i.e. no accretion, hence the con- 
stant baryon mass). Obviously the spin-down evolution 
is a function of time, therefore the x-axis of Figs. [H]could 
just as well be replaced by time. How the frequency (or 
equivalently the central density) varies with time will de- 
pend on the spin-down rate of the star. The computation 
of this quantity is not trivial, and depends on properties 
like magnetic field and/or gravitational radiation emis- 
sion. Those issues are beyond the scope of the current 
paper, and therefore we show our results as a function 
of frequency and central density, always keeping in mind 
the implicit time dependence. 

As shown in Figs. [S] and [51 when rotating at their Ke- 
pler frequency, neutron stars can attain higher masses. 
Loosely speaking this may be explained by the centrifu- 
gal force that provides an extra support against gravi- 
tational collapse. Even though we can find stable rotat- 
ing neutron stars with masses up to ^ 2.5 solar masses, 
these objects will collapse into black holes during spin- 
down evolution, which might be observable by the sudden 
stop of the neutrino signal originating from the star. As 
shown in Fig. [SJ there is no stable static neutron star with 
baryon mass greater than 2.18 solar masses predicted by 
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our model. 

The rotation also alters the redshift of the stars sig- 
nificantly. In Fig. [7] we show the forward redshift as 
a function of frequency for the stars of constant baryon 
mass of Fig. \5\ As one can see, the forward redshift is 
substantially modified as the star's rotational frequency 
is reduced. In the extreme case of high frequencies the 
forward redshift becomes negative. For comparison pur- 
poses we have also plotted the redshift of the equivalent 
stellar sequences with no quark phase (purely hadronic 
matter) . The redshift of these objects is almost identical 
to those of hybrid stars, with a very slight deviation for 
very high frequencies. 
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FIG. 7. (Color online) Forward redshift for the constant 
baryon mass sequences of Fig. [5] The dotted lines represent 
the redshift of the equivalent sequence for purely hadronic 
matter, which are almost identical to that of Hybrid stars. 

Just as the redshift is changed during the spin-down 
evolution, the moment of inertia of the star should be 
modified as well. The moment of inertia as a function of 
frequency for the stars shown in Fig. [5] are given in Fig. 

e 

Another particular feature of the model investigated 
here is the relatively low level of strangeness featured by 
the compact stars. As shown in Fig. |3]the A states (the 
first strange particle states to be populated) start to be 
occupied just before the onset of quark matter, which 
triggers the unstable branch of compact stars. Thus, the 
phase transition to quark matter suppresses the presence 
of strangeness in the object in the local charge neutral- 
ity case. Only for the densest objects we find a small 
region near the core with a low population of A states. 
The spin-down will also affect the strangeness of the ob- 
ject, increasing the radius of strangeness content near 
the core as the object spins down and becomes denser. 
This result is shown in Fig. [9] where the relative popu- 
lation {pA/pb) for the A states is shown for frequencies 
of 988.4 Hz (a), 746.0 Hz (b), 356.3 Hz (c) and 0.0 Hz 
(d). This result represents the strangeness content of the 



FIG. 8. (Color online) Moment of inertia for the constant 
Baryon mass sequences of Fig. [5] During the spin-down the 
geometry of the objects change significantly, and therefore the 
moment of inertia for these stars is substantially modified. 



star, since no other particles containing strangeness are 
present. Fig. [9] shows that for higher frequencies (and 
thus lower densities) strangeness states are either absent, 
or are very lowly populated. As the object spins down, 
and its density increases, the strange core expands and 
becomes more highly populated. For comparison pur- 
poses in Fig. [TU] we also show the strangeness content 
for a Mq/Msuu = 2.38 stars in the pure hadronic model. 
In this case, since there is no phase transition to quark 
matter, the stars are able to attain higher masses (as can 
be seen in Fig. [5]). Furthermore the hyperons are not 
suppressed by the onset of the quark phase, and these 
stars attain a higher strangeness content. As shown in 
Fig. [ini the A fraction at the core of the non-rotating 
star is 3.4 times higher than that found at the core of 
densest (stable) hybrid star. Once more we recall that 
the different frequencies shown in Figs. [H] and [TU] repre- 
sent different stages of time along the spin evolution of 
the star, with the time scale depending on the spin-down 
rate of the object. 



C. Cooling Process 

Another method of probing the inner core of compact 
stars is by investigating its thermal evolution. All the 
thermal processes taking place in a compact star strongly 
depend on its composition, therefore by comparing the- 
oretical predictions with observed thermal data, one can 
obtain valuable information regarding the core of neutron 
stars. 

The thermal evolution of neutron stars is dominated by 
neutrino emissions for the first 1000 years (maybe more 
in the slow cooling scenario) [l^, later being replaced 
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FIG. 9. (Color online) Strangeness content for a compact star with Mo/Msun = 2.18, for frequencies of 988.4 Hz (a), 746.0 Hz 
(b), 356.3 Hz (c) and 0.0 Hz (d). The z axis represents the relative population of A (pA/pb)- 



by photon emission from the surface. The direct Urea 
process [29] (DU process hereafter) is the most efficient 
coohng mechanism in a neutron star. With emissivities 
of the order of lO^^erg cm^ s~^ [l^, neutron stars in 
which the DU process takes place will cool very quickly. 
Due to momentum conservation however, the DU process 
can only take place when the proton fraction reaches a 
certain value (which depends on the underlying equation 
of state, but it is usually between 11 - 15 %). Compact 
stars, whose proton fraction are below this threshold, will 
feature a slower cooling. Other processes, like pairing and 
meson condensates, for example, might have an influence 
on whether or not a neutron star will feature a fast or slow 
cooling. These topics are beyond the scope of the current 
research and will be pursued in future investigations. 

We solve the thermal evolution equations ^] for the 
spherically symmetric stars shown in Fig. [5) The emis- 
sion processes considered for the core are the direct Urea 
process [l^, the modified Urea process [3l|, and the 



Bremsstrahlung process (sil . [s^ ; as for the crust we con- 
sider the electron Bremsstrahlung process ^] , electron- 
positron annihilation [s^ and plasmon decay [S^l . The 
results are shown in Fig. [11] where we show the cool- 
ing curves for neutron stars with different gravitational 
masses. The symbol Too denotes the redshifted temper- 
ature as detected by an observer at an infinite distance 
from the star. 

We have also plotted some pro minent observed tem- 
peratures of neutron stars (26l. l27j. Since there are two 
estimates for the age of pulsars, we have plotted two sets 
of observed data. The first set is for objects whose age 
estimate is based on the observed spin down rate, these 
are represented by squares. The second set is for ages 
obtained by tracking the pulsar back to its original su- 
pernova remnant (kinematic age, represented by circles). 

As shown in Fig. [TT] our model is in relatively good 
agreement with the observed data, with the exception 
of a few high temperature pulsars. It is possible that 





FIG. 10. (Color online) Strangeness content for a compact star with Mo/Msun = 2.38 in the pure hadronic model. The 
frequencies are 828.3 Hz (a), 639.9 Hz (b), 394.8 Hz (c) and 0.0 Hz (d). The z axis represents the relative population of A and 
S" ((pa + PE-/Pb)- 



these objects feature non-standard processes like pairing 
or some re-heating mechanism which would explain their 
high temperature. The shaded area represents all possi- 
ble cooling curves for our model, that can be obtained 
by solving the cooling equations for the whole spectrum 
of stable hybrid stars shown in Fig. [5] We can also see 
that within our model, neutron stars with masses higher 
than 1.22 solar masses, which feature an enhanced cool- 
ing, cool down too fast to be in agreement with observed 
data. It is important to mention that these results can 
be further improved if one consider more sophisticated 
processes like pairing and meson condensation. 

Finally we have also investigated how spin down effects 
might affect the thermal evolution. As shown in Fig. 
[9l the increase in density that follows spin down has a 
strong effect on the strangeness content of the star. It 
is only natural to expect that other particle states will 
also be altered by the spin down. Particularly interesting 



is the proton, electron and neutron populations, since 
these particle population will dictate whether or not the 
direct Urea process takes place. This process can only 
take place if the following triangle inequality (and cyclic 
permutations of it) is satisfied 

ki+kf> ki, (8) 

with ki denoting the Fermi momentum of particle i. We 
calculated the radius threshold for different frequency 
stars, from the sequence with il/g /A/gun = 1.6. This re- 
sult is shown in Fig. I12[ where the gray shaded repre- 
sents the region in which the DU process is allowed to 
take place. 

The results of Fig. [T^ show that the region in which 
the DU process takes place changes substantially during 
the evolution of the compact star. For the case studied, 
for higher frequencies (where the density is lower) the DU 
process can only take place in a very small region near 
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IV. CONCLUSION 
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FIG. 11. (Color online) Cooling of a range of spherically sym- 
metric neutron stars from Fig. [21 The symbol Too denotes 
the redshifted temperature as detected by an observer at an 
infinite distance from the star. The onset of the direct Urea 
process happens for neutron stars with 1.22 solar masses. Any 
star with a lower mass will have a slow cooling. For masses 
above this value an enhanced cooling is achieved. Also shown 
is some prominent observed temperatures of neutron stars 
[2^.[27l|. Squares denote that the age was estimated based on 
the spin down rate, and circles indicate age estimates based 
on the motion of the pulsar with respect to its originating 
supernova remnant. 



the core. This regions grows as the stellar frequency de- 
creases and the object becomes denser. This hints that 
the coohng of compact stars might be slower if one con- 
sider effects of spin down. Current efforts are in progress 
to better understand this issue, and wih be discussed in 
a future paper. 

Finally we have also generated a diagram il{M) that 
shows the domain of frequencies and mass for which the 
direct Urea process is allowed. This result is shown in 
figure I13[ where any star with mass and frequency that 
falls within the dark gray shaded area allows for the di- 
rect Urea process to take place. The light grey shaded 
area, on the top of the diagram, represents a forbidden 
region, where no stars can be found since their frequency 
would be above their Kepler frequency. The white area 
represents stars in which the direct Urea process can- 
not take place. Currently the spin of many X-ray burst 
sources is known, and in some cases their masses can also 
be inferred (see (sB - fSTj and references therein). In the 
case of transient neutron stars, the core temperature can 
be inferred from the quiescent emission state [38]. Fur- 
thermore, for some neutron stars the spin frequency and 
the core temperature are bound (see for instance [39|). 
In the event that one might estimate the core tempera- 
ture, spin frequency and stellar mass, the diagram shown 
in Fig. 1131 might be used as a further test for this model. 



A major advantage of our work, when compared to 
other studies of hybrid stars, is that because we have only 
one equation of state for different degrees of freedom we 
can study in detail the way in which chiral symmetry is 
restored and the way deconfinement occurs in the stars. 
Such phenomena happens, for example, during the star 
spin down. 

We have found that SU(3) non-linear sigma model is 
suitable for the description of hybrid stars. The predicted 
maximum masses and the respective radii lie in the ob- 
served range. For a static object, our model predicts 
a star containing 2 km of hybrid matter (radius) sur- 
rounded by hadronic matter. In the event that the object 
is rotating at its Kepler frequency, the hybrid core be- 
comes an oblate ellipsoid with equatorial and polar radii 
of 1.18 and 0.87 km, respectively. The reduction of the 
quark core for a rotating object should not be surpris- 
ing, if one considers the reduction in density that follows 
rapid rotation. 

We have also investigated the cooling of spherically 
symmetric (static) neutron stars, whose composition is 
described by our model. We have found that the thresh- 
old for the direct Urea process is reached for stars with 
mass greater than 1.0 M©. We have compared the cool- 
ing curves predicted by our model with some prominent 
observed compact stars temperature. We have found 
that, within our model, objects with masses up to 1.22 
Mq are in good agreement with the observed data. Any 
star with mass above this value features a thermal evolu- 
tion that is too fast to be in agreement with the observed 
data. This result might seem inconsistent with the ob- 
served data, since most of these objects are expected to 
have a higher mass than 1.22 Mq. It is important to 
notice that the results shown in Fig. [TT] were obtained 
assuming a "froze- in" structure/composition. As noted 
in reference [38^, many observed accreting neutron stars 
are in agreement with the slow cooling scenario. The 
fact that they are accreting implies that their structure is 
changing, as their frequency is modified as a result of the 
accretion process. As we showed in Fig. [T^l the thresh- 
old radius for the direct Urea process strongly depends 
on the frequency of the star. Therefore, the fast cooling 
for stars with masses above 1.22 Mq might be deceiving, 
since if we consider rotation we might obtain stars with 
Mq = 1.6Mq without presenting the direct Urea process 
(if the frequency is high enough) . We could also see that 
a few observed objects present a very high temperature, 
that cannot be explained by our model. Most likely there 
is some non-trivial heating process taking place in these 
objects, which would explain why they are so warm at 
their relatively old ages, this topic however is beyond the 
scope of this work. It is important to mention that the 
quarks found in the 2 km hybrid core, have little effect on 
the cooling of these objects. Not only is the hybrid core 
very small, but the quarks are present at a smaller ratio 
than the hadrons, which allows the hadrons to dominate 
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FIG. 12. (Color online) Direct Urea threshold for a compact star with Mo/Msun ~ 1.6, for frequencies of 792.4 Hz (a), 656.8 
Hz (b), 496.5 Hz (c) and 0.0 Hz (d). The grey shaded areas denotes the region in which the direct Urea process is allowed to 
take place (denoted by z=l). 



the cooling processes. On the other hand the suppression 
of hyperons, due to the onset of the quark phase is impor- 
tant. Since the hyperons do not appear in great quanti- 
ties, some coohng channels (Hyperon Urea processes for 
instance) are not open. These channels are not very effi- 
cient cooling mechanisms [i^, nonetheless their absence 



slow down the cooling. 
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